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Abstract 

We calculate the 0(a s ) virtual corrections to the decay width for b — > sg 
in the standard model (g denotes a gluon). Also the corresponding 0{a s ) 
bremsstrahlung corrections to b — > sg are systematically calculated in this 
paper. The combined result is free of infrared and collinear singularities, in 
accordance with the KLN theorem. Taking into account the existing next-to- 
leading logarithmic (NLL) result for the Wilson coefficient C| ff , a complete 
NLL result for the branching ratio i3 NLL (6 — > sg) is derived. Numerically, we 
obtain ,g NLL = (5.0 ±1.0) x 10 -3 , which is more than a factor of two larger 
than the leading logarithmic result £> LL = (2.2 ± 0.8) x 10 -3 . 
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I. INTRODUCTION 



The theoretical predictions for inclusive decay rates of 5-mesons rest on solid gounds 
due to the fact that these rates can be sytematically expanded in powers of (Aq CD /m b ) 
where the leading term corresponds to the decay width of the underlying 6-quark decay As 
the power corrections start at 0(AQ CD /m(,) 2 only, they affect these rates by at most 5%. 
Thus the accuracy of the theoretical predictions is mainly controlled by our knowledge of 
the perturbative corrections to the free quark decay. 

The charmless inclusive decays, B — > Xj, where X^ denotes any hadronic charmless final 
state, are an interesting subclass of the decays mentioned above. At the quark level, there 
are decay modes with three-body final states, viz. b — > q'q'q, (q' — u, d, s; q — d, s) and the 
modes b —>■ qg, with two-body final state topology, which contribute to the charmless decay 
width at leading logarithmic (LL) accuracy. Next-to-leading logarithmic (NLL) corrections 
to the three-body decay modes were started already some time ago in ref. ||, where radiative 
corrections to the current-current diagrams of the operators 0\ and O2 were calculated, 
together with NLL corrections to the Wilson coefficients. Later, Lenz et al. included the 
contributions of the penguin diagrams associated with the four-Fermi operators 0\, ■■■,Oq] 
the effects of the chromomagnetic operator 0§ were taken into account to the relevant 
precision needed for a NLL calculation. Up to contributions from current- current type 
corrections to the penguin operators, the NLL calculation for the three quark final states is 
complete. 

In the numerical evaluations of the charmless decay rate, the two body decay modes 
b — > qg were added in refs. [|,|J at the LL precision, as the full NLL predictions were 
missing. It is exactly this gap which we try to fill in the present letter. We will present 
the results of a calculation for the branching ratio B(b — > sg) where NLL corrections are 
systematically included. This implies that besides virtual corrections to b — > sg also the 
process b — > sgg has to be taken into account, as it gives contributions at the same order 
in perturbation theory. The LL prediction for the branching for b ^ sg is known to be 
B(b — ► sg) ~ 0.2% M; also the process b — ► sgg has been studied in the literature J7|||. 
In p[ a complete calculation was performed in regions of the phase space which are free of 
collinear an infrared singularities, leading to a branching ratio for b — > sgg of the order of 
10~ 3 . A complete NLL calculation requires, however, a regularized version for the decay 
width F(b — > sgg) in which infrared and collinear singularities become manifest. Only 
after adding the virtually corrected decay width T(b — > sg) a meaningful physical result is 
obtained. In addition, as we will see later, also the tree-level contribution of the operator 
Os to the decays b — > sff, with f — u,d, s, has to be included. 

The decay b — > sg gained a lot of attention in the last years. For a long time the 
theoretical predictions for both, the inclusive semileptonic branching ratio and the charm 
multiplicity n c in 5-meson decays were considerably higher than the experimental values . 
An attractive hypothesis, which would move the theoretical predicitions for both observables 
into the direction favored by the experiments, assumed the rare decay mode b — > sg to be 
enhanced by new physics. 

After the inclusion of the complete NLL corrections to the decay modes b — > cuq and 
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b — > ccq (q = d, s) [p70| , both the CLEO- and the LEP-data |TT[ are now in agreement with 



theory |To| , [T2[| , if one allows the renormalization scale \i to become as low as mb/A. We would 
like to stress that there is still some room for enhanced b — > sg, in particular when using 
higher values for the renormalization scale. For theoretical motivations of enhanced b — > sg, 



see ref. [|13[] and references therein. 

We also would like to mention that the component b — > sg of the charmless hadronic 
decays is expected to manifest itself in kaons with high momenta (of order m&/2), due to 



its two body nature Some indications into this direction were reported by the SLD 

collaboration |T5|| . For overviews on enhanced b — ► sg, see e.g. refs. |TB| , PT| . 

The remainder of this letter is organized as follows: In section [nj we briefly review the 
theoretical framework. Section [TTI| is devoted to the calculation of the virtual corrections to 
the decay width for b —>■ sg, while section deals with the calculation of the bremsstrahlung 
corrections to b — > sg. In the short section [TV] in between, the decay width for the tree level 
processes b — > sff mentioned above, is given. In section [Vj we give the expressions for the 
NLL branching ratio i3 NLL (5 — > sg), which combines the processes b — > sg, b — > sgg and 
6 — > s//. Finally, in section |VII| the numerical results for £> NLL (6 — > s#) are presented. 



II. THEORETICAL FRAMEWORK 

The analysis of the decays b ^ sg and 6 — > sgg starts with introducing the effective Hamil- 
tonian 

K ff = --^aj;^), (i) 

where V^j are elements of the CKM matrix, Oj(/i) are the relevent operators and Cj(//) are 
the corresponding Wilson coefficients. The full set of operators needed for our application, 
can be seen in ref. ||18|1 . As the Wilson coefficients of the gluonic penguin operators are 
small (see eq. (30) in ref. fl8|), we neglect them when calculating radiative corrections; we 



therefore only list the explicit form of the operators 0\, O2 and Os : 

Oi = (s L ^T A c L ) (c L ^T A b L ) , 2 = (sl^cl) (c L rh) , 

08 = ^ m 6 (//) (s L a^T A b R ) G% . (2) 

Here T A stand for the 5?7(3) co i or generators. The small CKM matrix element V u b as well as 
the s-quark mass are also neglected. 

It is well-known that in this formalism the large QCD logarithms, present in the decay 
amplitude for b — > sg, are contained in resummed form in the Wilson coefficients when 
choosing the renormalization scale \i at the order of m&. The LL (NLL) Wilson coefficients 
contain all terms of the form a™\n n (mb/M) (a s a™\n n (mb/M)), where M = m t or raw an d 
n = 0,1,2,.... 

The LL prediction for the decay amplitude for b — > sg is then obtained by calculating 
the matrix elements (sg\Oi\b) at order g s and weighting them with the leading logarithmic 
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FIG. 1. Complete List of two-loop Feynman diagrams for b — > sg associated with the operators 
0\ and 02- The fermions (b, s and c quarks) are represented by solid lines; the wavy lines represent 
gluons. The crosses denote the possible locations where the gluon is emitted. 



Wilson coefficients. In the naive dimensional regularization (NDR) scheme which we use in 
this paper, there are one-loop contributions of order g s for % = 3,4,5,6 and the tree-level 
contribution of Og. The effect of the matrix elements for i = 3, 4, 5, 6 can be absorbed into 
the effective Wilson coefficient (see 0) Cf = C$ + C 3 - | C 4 + 20 C 5 - f C 6 . 

The NLL corrections for the decay amplitude for b —> sg receives two contributions: 
The first one arises when combining the lowest order matrix elements (of order g s ) with 
the NLL Wilson coefficient, while the second one arises when calculating explicit order a s 
corrections to the matrix elements of the operators which are then weighted with LL Wilson 
coefficients. As the operators 0\ and O2 have vanishing matrix elements for b — > sg at 
order g s and the NLL corrections connected with the operators O3, Oq are neglected, the 
only Wilson coefficient needed to NLL precision is that of the operator Og. The necessary 
ingredients, the anomalous dimension matrix to 0(a 2 s ) and the 0(a s ) matching condition for 
the operator Og were given in refs. and (|TjJ , respectively. A practical parametrization for 
Cg ff (/i) will be given in ref. pOfl . A rather complete list of references on matching conditions, 
anomalous dimension matrices, and on the process b — > 57, which is similar in many respects 



to b — ► sg, is given e.g. in ref. [21 



III. VIRTUAL CORRECTIONS TO O x , 2 AND O, 



A. Virtual corrections to 0\ and O2 



As the one-loop matrix elements of the operators 0\ and O2 vanish, we immediately 
turn to the two-loop contributions. A complete list of Feynman diagrams for the matrix 
elements (sg\Oi\b) (i = 1,2) is shown in fig. |I| The diagrams in fig. [T] a), b), c) and d), 
in which the emitted gluon is replaced by a photon are the relevant diagrams for b — > 57; 
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these were calculated in ref. [22|. The three possible diagrams in fig. [1] g) cancel each other 



when considering the process b — *> S7; however, they give a non-vanishing contribution^ to 
b — > sg. There is also another difference: while each of the figures a), b), c) and d) forms a 
gauge invariant subset in b — > sj, this is no longer true for b — * sg; a gauge invariant result 
is only obtained when all the diagrams in fig. [I] are summed. 

The various two-loop integrals are calculated by the standard Feynman parameter tech- 
nique. The heart of our procedure, which is explained in detail in ref. |22[ for one of the 
diagrams in fig. [3] a), is to represent the rather complicated denominators in the Feynman pa- 
rameter integrals as complex Mellin-Barnes integrals |23|| . After inserting this representation 
and interchanging the order of integration, the Feynman parameter integrals are reduced 
to well-known Euler Beta-functions. Finally, the residue theorem allows to represent the 
remaining complex integal as the sum over the residues taken at the pole positions of Beta- 
and Gamma-functions; for a generic diagram G, these steps naturally lead to an expansion 
in the ratio z = (m c /mb) 2 of the form 



2 

m 



G = c + J2cnmZ n L m ; z = -§ ; L = \nz, (3) 



n,m 



ni 



where the coefficients Cq and c nm are independent of z. The power n in eq. (|^) is in general 
a natural multiple of 1/2 and m is a natural number including 0. In the explicit calculation, 
the lowest n turns out to be n = 1. This implies the important fact that the limit m c — > 
exists. As was shown in |[22| , the power m of the logarithm is bounded by 4, independently 
of the value of n. In our results, which we will present below, all terms up to n = 3 are 
retained. 

We first present the final result for the dimensionally regularized matrix element M 2 = 
(sg\02\b) which represents the sum of all two- loop diagrams in fig. Ill: 




(4) 



M 2 = -^(sg\0 8 \b) tTee 

47T 

where the real- and imaginary part of r 2 read 

Re(r 2 ) = ^{ - 2170 - 54tt 2 + z[48816 - 252tt 2 + (22680 - 1620tt 2 )L 
+2808L 2 + 612L 3 - 6480C(3)] 

-12672ttV /2 + 2 2 [66339 + 1872tt 2 + (-40446 + 1512tt 2 )L 
+6642L 2 - 1008L 3 + 7776C(3)] 
+z 3 [-3420 - 60tt 2 - 6456L + 7884L 2 ]} 

Im(r 2 ) = — { - 28 + ^[549 - 24tt 2 + 153L + 72L 2 ] 

+z 2 [-432 + 30tt 2 + 54L - 90L 2 ] + z 3 [-259 + 192L] } . (5) 



We thank M. Neubert for making us aware of these diagrams. 
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The symbol ( stands for the Riemann Zeta function, with £(3) ~ 1.2021. Finally, the 
quantity (sg\0 8 \b) tTee denotes the tree level matrix element of the operator 8 . 

To obtain the renormalized matrix element M 2 en associated with the operator 2 , the 
corresponding counterterms have to be included. This means that we have to take into 
account the one-loop matrix elements of the four Fermi operators 5Z 2 jOj (j = 1, ...,6) and 
the tree level contribution of the magnetic operator 5Z 28 8 . In the NDR scheme the only 
non- vanishing contributions to b —>■ sg come from j = 4, 8. The operator renormalization 
constants Zij can be extracted from the literature |Tj| in the context of the leading order 



anomalous dimension matrix. One obtains the counterterm contribution 

Mf = {sg\5Z 2A A + SZ 28 O s \b) = (-£- - ( ^ + - ^ ~ ] (sg\0 8 \b) tiee . (6) 

\ 3o7r e \ ji j 7r 108 e / 

We note that there are no one-loop contributions to the matrix element for b —>■ sg from 



counterterms proportional to the evancesent operator P 12 given in appendix A of ref. |18 
Adding the regularized two-loop result in eq. (f|) and the counterterm in eq. (^|), we find 
the renormalized result for M 2 in the NDR scheme: 

Mr = (sg\0 8 \b) tice g [i 2 In ^ + r 2 ) , (7) 

where r 2 is given in eq. (|5|) and £ 2 = 70/27. 

By doing analogous steps, we obtain the renormalized version of Mi = (sg\Oi\b): 

Mr = (sg\0 8 \b) tiee g (li In ^ + ri \ , (8) 

with l x = 173/162 and 

Re(ri) = — {4877 - 54tt 2 + 36z[1086 + 29vr 2 + (360 + 36tt 2 )L 

3888 

+51L 2 + 8L 3 + 144C(3)] 

-12672ttV /2 + 9^ 2 [6615 - 80tt 2 + (-4494 + 384vr 2 )L 
+864L 2 - 148L 3 + 864C(3)] 
+12z 3 [93 + 76tt 2 - 1186L + 900L 2 ]} 

Im(ri) = -^{25 + 6z[75 + tt 2 + 24L - 3L 2 ] 

+6^ 2 [-171 + 19tt 2 + 72L - 57L 2 ] + 2z 3 [-421 + 192L]} (9) 

For z > 1/4 the imaginary parts of r\ and r 2 must vanish exactly; our results fulfill this 
property to high accuracy when retaining terms up to z 3 in the expansion. 
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c) 
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d) 



08^ 
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e) 



08 i 



FIG. 2. Complete list of Feynman diagrams associated with the operator The analogues 
of a), c) and d), where the virtual gluon hits the 6-quark instead of the s-quark, are not shown 
explicitly. The real gluon can be attached to any of the crosses shown in a). 



B. Virtual corrections to Os 



We present now the results of the virtual corrections to the matrix element M§ = 
(sg\0$\b). As the contributing Feynman graphs in fig. ||| are all one loop diagrams, the 
computation of Mg is straightforward. We use dimensional regularization for both, the ul- 
traviolet and the infrared singularities. Singularities which appear in the situation where the 
virtual gluon becomes almost real and collinear with the emitted gluon are also regulated 
dimensionally; on the other hand, those singularities where the almost real internal gluon 
is collinear with the s-quark, are regulated with a small strange quark mass m s ; the latter 
manifest themselves in logarithmic terms of the form ln(p), where p = {m s jm^) 2 . 

We were able to separate the ultraviolet 1/e poles from those which are of infrared 
(and/or collinear) origin. For ultraviolet poles we use the symbol 1/e in the following, while 
collinear and infrared poles are denoted by 1/erR. 

As the results of the individual diagrams are not very instructive, we only give their sum: 



47T 



f 8 (sg\0 8 \b) 



tree j 



(10) 



with 



fs 



(41n(p) + 9 + 9in) 



11 

3e" 



1 

3 



59tt 5 
~Y1 



+ 1 



3em 

-81n(p) + 21n 2 (p) 




Mr 



-2c 



+ 



(11) 



An ultraviolet finite result is obtained by adding the contribution from the counterterm 
which is generated by expressing the bare quantities in the tree-level matrix element of 0§ 
by their renormalized counterparts. It has the structure 



M* = 6R(sg\0 8 \b} tme 



(12) 



where the factor 5R is given by 5R = J Z 2 (rrib) \/Z : 



Z% Zg s Z 8S — 1 . 
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Ziiyn-b), Z2(m s ) and Z3 denote the on-shell wave function renormalization factors of the 
6-quark, the s-quark and the gluon, respectively. Z ga and Z mb denote the MS renormalization 
constants for the strong coupling constant g s and the 6-quark mass, which appear explicitly 
in the definition of the operators (see eq. (0)). Finally, Z 88 is the renormalization factor 
of the operator Os- The Z-factors of the fields, of the masses and of the strong coupling 
constant are given in text books, while Z§g can be extracted from the anomalous dimension 
matrix in ref. |Tj|; we therefore immediately give the expression for 5R: 



5R=- 



An 



11 

3e 6ei R 



31 _ 81n !^_2 Eln 



VTit 16 
jj, 3 



2 hip 



(13) 



The term / originates from fermion self-energy diagrams contributing to the on-shell renor- 
malization constant Z% of the gluon field; / runs over the five flavors u, d, s, c and b. 

When adding the regularized matrix element of 0§ in eq. (|ToD and the counterterm 
contribution Mf* in eq. (|i~2"|), we obtain the renormalized result 



M 8 ren 



An 



fr(sg\o 8 \b) 



tree ) 



(14) 



with 



,-r 
J8 



(8 ln(p) + 49 + 18i7r) 



-IR 



6em 



59tt 2 




-2c 



29, m b 
— m — 

3 p 



- In p H — In p m . 

3^3 ^3 



(15) 



We anticipate that the singular terms of the form l/ej R , 1/erR and lnp in eq. (|T^) will 
cancel against the corresponding singularities in the result for the gluon bremsstrahlung 
corrections to b — >• sg. On the other hand, the terms ln(m///i), which also represent some 
kind of singularities for the light flavor f — u,d,s are not cancelled in this way. Keeping 
in mind that they originate from the fermionic contribution to the renormalization factor 
Z 3 , it is expected that they will cancel against the logarithms present in the decay rate 

T(b^8ff)(f=U,d,8). 



C. Virtual corrections to the decay width T(b — > sg) 



We are now in the position to write down the renormalized version of the matrix 
M ren (b — > sg) element for b — ► sg, where the virtual order a s corrections are included. 
We obtain: 



M ren (6 -> sg) 



AG pi 



v; s v tb 



cf + ^ 



a. 
An 



Cl{h In ^ + n) + C° 2 (£ 2 In ^ + r 2 ) + 



/^0,eff rr 
8 J8 



(sg\Os{n)\b) 



tree • 



(16) 
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The quantities n, r 2 and /g en are given in eqs. @, @ and (|T5|), respectively. As eq. (|16D 
shows, Cg ff is the only Wilson coefficient needed to NLL precision. For the following it is 
useful to decompose it as 



8-8 4vr 8 



(17) 



The symbol {sg\Os(^)\b) tree in eq. (|16|) denotes the tree level matrix element of 8 (/i), which 
contains the running 6-quark mass and the strong running coupling constant at the scale fi. 
In order to get expressions where the 6-quark mass enters as the pole mass, and the strong 
coupling constant enters as g s (m b ), we rewrite (sg\Os(fi)\b) tr ee as 



(sg\0 8 (n)\b) tree = {sg\0 8 \b) t 



2a s m b 4a s a s m b 
1 + In -— + — A, hi — 

7T jl 6 TX 47T jJL 



A> = f ■ ^ 



The symbol (sg\0$\b) tTee then stands for the tree level matrix element of 8 in which m b (ii) 
and g s have to to be indentified with the pole mass m b and g s (m b ), respectively. After 
inserting eqs. ( |17|) and (p~8|) into eq. (0), the corresponding decay width r virt is obtained 
in the standard way. One gets: 



24vr 4 



\GFV t * s V tb \ 2 



2 a s n 0,eS 



47T 



2C 8 1 ' eff + 2(8 + /3 )ln^C; 



^C 8 °' eff + 2 C?& In ^ + He(n)) + 2 C 2 °(£ 2 In ^ + Re(r 2 )) 



+2 C 8 °' eff Re(/ 8 ren ) 



-2e 



l + 2e 



- 16 )e >)] } 



(19) 



We note that due to the infrared poles present in /g en the phase space integrations were 
done consistently in d = 4 — 2e dimensions. 



IV. Og CONTRIBUTION TO THE DECAY WIDTH T(b s//) 



As discussed at the end of section |III B| , we should take into account the contribution 
of the operator 8 to the process b — > sf f (/ = u, d, s), in order to cancel the unphysical 
logarithms of the form ln(m///z) in the virtual corrections to b — > sg. The 8 contribution 
to the decay width T 8 (b — > sf f) yields 



Ts(b - sff) 



^ ml\G F V*V tb C° 8 ' eSl2 



72tt 5 



In 



mi, 
2m/ 



(20) 



V. GLUON BREMSSTRAHLUNG CONTRIBUTIONS 



In this section we discuss the gluon bremsstrahlung corrections to b — > sg, i.e. the process 
b — > sgg. As in the case of the virtual corrections, we neglect contributions from the gluonic 
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b 08 s b 08 s b 08 s 

(8 — ■- — & 



^ 1° °?r~ 




5' e. 



FIG. 3. Complete list of Feynman diagrams for b — > 555, associated with the operators Oi, O2 
and 8 . 

penguin operators as their Wilson coefficients are rather small. In this approximation, the 
matrix element M brems (6 — > sgg) is of the form 



AG F i 



Cl M brems + C 2 ° M 2 brcms + C 8 °' eff M brems , (21) 



where the three terms on the r.h.s. correspond to the contributions of the operators 0±, O2 
and 8 , respectively. The corresponding Feynman diagrams are shown in fig. [3[ We note 
that in eq. Q2"T| ) only the leading order pieces of the Wilson coefficients are needed. 

The decay width r brems (6 — > sgg) is then obtained by squaring M brems , followed by 
phase space integrations. These integrals are plagued with infrared and collinear singulari- 
ties. Configurations with one gluon flying collinear to the s-quark are regulated by a small 
strange quark mass m s , while configurations with two collinear gluons, or one soft gluon 
are dimensionally regularized. As in the calculations of the virtual corrections, we write the 
dimension as d = 4 — 2e. (Note that e has to be negative in order to regulate the phase 
space integrals). 

When squaring M brems in eq. ([H]), nine terms are generated, which we denote for 
obvious reasons by (O^O*), (0^0*), (O^O*), (0 2 ,0*), (0 8 ,0*). We find that all 
these quantities are free of infrared and collinear singularities, except (0 8 ,0|). Hence, one 
can put m s = in the finite terms and evaluate the phase space integrals in d = 4 dimensions. 
In the following, we denote this finite contribution to the decay width by rg^ cms . It turns 
out that only ~ 5% of the total NLL correction are due to rg^ cms . As the analytical results 
for this finite piece, written in terms of two-dimensional integrals, are rather lengthy, we 
skip the explicit expressions in this letter; we stress, however, that Fg^ ems , although small, 
will be retained in the numerical evaluations. 

We now turn to the (0 8 , 8 ) contribution, denoted by r 88 ems . After a lengthy, but 
straightforward calculation, we obtain (p = (m s /mj,) 2 ; V = a s ml \G FV t *V t b\ 2 / (24:7T 4 )) 



1 ss — 



-it 

"18 67 + 8 hip l2 , A 62vr 2 

— + - 4 In 2 p + 12 In p + 240 



22) 



12 71 L e m e m 

The total decay with for b — > sgg is then given by 

r brcms (6 -> sgg) = r b ^ ems + r b ^ ems . (23) 
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VI. COMBINED NLL BRANCHING RATIO FOR b -» sg AND b -» 



In this section we combine the decay widths for the virtually corrected process b — > sg 
and the bremsstrahlung process b — > sgg to the total NLL decay width decay T NLL (6 — > sg). 
We also absorbe in this quantity the 0$ induced contribution to the processes b — > sf f 
(/ = u, d, s), as discussed in section fV[ From the explicit formulas for r virt and r brems one 
can see that the infrared and collinear singularities cancel in the sum. The terms containing 
logarithms of the light quark masses rrif, present in the result for r virt , are cancelled when 
combined with T 8 (b — > sff). Putting together the individual pieces, the final result for 
r N (b — > sg) can be written as 



^NLL 



(b -> sg) 



24tt 4 



G F V*V tb \ 2 \D\ 2 + T 



-ibrcms 
fin 



(24) 



with 



Air 



a 



l.eff 



16 



— C 8 ' + Ciyii log h n 



+C° 2 [£ 2 log ^ + r 2 ] + C 8 °' eff [(4 + 8 + Po) log ^ + 



(25) 



A remark concerning the modulus square of the function D is in order: By construction, 
this square is understood to be taken in the same way as the in the virtual contributions, i.e. 
by systematically discarding the 0(a 2 s ) term. In this sense, the quantity D can be viewed 
as an effective matrix element. We stress however that, besides the virtual corrections, also 
the informations of rgg ems and T s (b — > //) are contained in the function D. 

The quantities r x and r 2 appearing in eq. (|25|) are given in eqs. (||) and ©, respectively. 
The explicit expressions for £i, £ 2 , is and r 8 , read 



173 
162 



70 
27 



19 
"3 



1 

18 



351 - 19vr z 



777 

36 ln2 + 61n^ 

ml 



(26) 



Note that all scale dependent quantities in eqs. (p4|) and (p5|) are understood to be evaluated 
at the scale /x, unless indicated explicitly in the notation. 



2 and (is + 8 + /?o) are identical to the anomalous 

0,cff 



We would like to point out that i\ 
dimension matrix elements 7i 8 e£f , 728 eff ? an d ls^^ respectively, which are given in ref. |I8"| 
This is what has to happen: Only in this case the leading scale dependence of Cg' eff (/i) gets 
compensated by the second term in eq. (|25|). 

The NNL branching ratio B NLL (b — * sg) is then obtained as 



B (b -> sg) 



»exp 



(27) 



where B°i P denotes the experimental semileptonic branching ratio of the I?-meson. T s i stands 
for the theoretical expression of the semileptonic decay width of the 5-meson. Neglecting 
non-perturbative corrections of the order (AQ CD /mf,) 2 , T s \ reads (with x c = (m c /m b )) 
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FIG. 4. Scale (/i) dependence of the function D (see eq. (p5|)) in various approximations: The 
long-dashed line shows Cg' e ; the short-dashed line corresponds to putting r± = r 2 = r§ = 0; the 
dotted line is obtained by only putting r 2 = 0; the solid line shows the full function D. See text. 



_ . G 2 F m 5 
r s i «1(M cev e ) - 
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l + ^-h sl (x c ) + 0(a 2 s ) 

Z7T 



(28) 



with the phase space function g(x c ) = 1 — 8x^ — 24 x^ lnx c + — x®. The analytic 
expression for h s \(x c ) can be found in ref. [p4 ]. The approximation (taken from ref. Q) 



h sl ( Xc ) = -3.341 + 4.05 (x c - 0.3) - 4.3 (x c - 0.3)' 



(29) 



which we use in the following, holds to an accuracy of one permille for 0.2 < x c < 0.4. 



VII. NUMERICAL RESULTS FOR THE COMBINED BRANCHING RATIO 



We first discuss the sizes of the various NLL corrections at the level of the function 
D, defined in eq. (p5|). As already stated, the terms containing explicit logarithms of the 
ratio {rrib/ n) get compensated by the scale dependence of the first term on the r.h.s. of eq. 
(p5|). This feature can be observed in fig. ||, when comparing the two dashed lines. The 
long-dashed line represents only the first term C§ of the function D, while the short-dashed 
line shows D, in which r 1; r 2 and r 8 are put to zero. As expected, the short-dashed line 
has a milder /i-dependence. When switching on also r\ and r 8 (but keeping r 2 = 0), the 
resulting curve, shown by the dotted line, stays close to the short-dashed curve and the scale 
dependence remains very mild. However, when switching on also r 2 , the situation changes 
drastically. The resulting solid line, which represents the full NLL D function, implies that 
the term containing the two- loop quantity r 2 , induces a large NLL correction. As this large 
correction term contains a factor a s {fi) C 2 (/z), it is of no surprise, that the NLL prediction 
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FIG. 5. Branching ratio B(b —* sg) as a function of the scale /i in various approximations: The 
dashed and the solid lines show the LL and the NLL predictions, respectively; the dotted line is 
obtained by putting t\ = r 2 = r% = r{j° ems = in the NLL expression for T NLL (6 — ► sg) in eq. 
(Mj. See text. 



for the function D suffers from a relatively large scale dependence, as illustrated by the solid 
line. 

The NLL branching ratio i3 NLL (6 — > sg) is then obtained as described in section [VI]. The 
result is shown by the solid line in fig. [|. For the input values, we choose: m b = (4.8 ± 0.2) 
GeV, im c /m b ) = (0.29 ± 0.02), a s {m z ) = 0.119 ± 0.003, \V*V tb /V cb \ 2 = 0.95 ± 0.03, B™ p = 
(10.49 ± 0.0046)%, and mf ole = (175 ± 5) GeV. As the scale dependence is rather large, we 
did not take into account the error due to the uncertainties in the input parameters. Based 
on fig. H|, we obtain the NLL branching ratio 

B NLL {b^sg) = (5.0 ±1.0) x 10~ 3 , (30) 

which is more than a factor two larger than the LL value 

B LL (b^ sg) = (2.2 ±0.8) x 10~ 3 , (31) 

extracted from the dashed line in fig. [5| As stressed in the discussion of the function D, 
the main enhancement is due to the virtual- and bremsstrahlung corrections to b —* sg, 
calculated in this paper. At the level of the branching ratio, this fact is illustrated by the 
dotted line in fig. |, which is obtained by discarding r^ ems and by switching off ri, r 2 and 
r 8 in the expression for r NLL (6 — > sg) (see eq. (p4"D). 

The largest uncertainty due to the physical input parameters (other than n) on B(b — > sg) 
results from the charm quark mass. Varying x c = m c /m b between 0.27 and 0.31 and choosing 
H = m b , the resulting uncertainty amounts to ~ ±6%. 
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